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Abstract. We show that aperiodic and quasiperiodic variability of bright LMXBs - atoll and Z- sources, on ~ sec - msec 
time scales is caused primarily by variations of the luminosity of the boundary layer. The emission of the accretion disk is 
less variable on these time scales and its power density spectrum follows Pdisk(/) oc f~^ law, contributing to observed flux 
variation at low frequencies and low energies only. The kHz QPOs have the same origin as variability at lower frequencies, 
i.e. independent of the nature of the ’’clock”, the actual luminosity modulation takes place on the neutron star surface. The 
boundary layer spectrum remains nearly constant in the course of the luminosity variations and is represented to certain 
accuracy by the Fourier frequency resolved spectrum. In the investigated range of M ~ (0.1 — l)MBdd it depends weakly 
on the global mass accretion rate and in the limit M ~ M^dd is close to Wien spectrum with kT ^2 .4 keV. Its independence 
on the global value of M lends support to the theoretical suggestion bv llnogamov & SunvaevI il999ll that the boundary layer 
is radiation pressure supported. 

Based on the knowledge of the boundary layer spectrum we attempt to relate the motion along the Z-track to changes of 
physically meaningful parameters. Our results suggest that the contribution of the boundary layer to the observed emission 
decreases along the Z-track from conventional ~ 50% on the horizontal branch to a rather small number on the normal 
branch. This decrease can be caused, for example, by obscuration of the boundary layer by the geometrically thickened 
accretion disk at M ~ MEdd. Alternatively, this can indicate significant change of the structure of the accretion flow at 
M Msdd and disappearance of the boundary layer as a distinct region of the significant energy release associated with the 
neutron star surface. 
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1. Introduction 

Accreting neutron stars in low mass X-ray binaries (LMXB) 
are among the most luminous compact X-ray sources in the 
Milky Way. A number of them have luminosities exceeding 
~ few X 10^® erg/s and presumably accrete matter at the level 
close to the critical Eddington accretion rate. In the bright 
state these sources have rather soft X-ray spectra, indicating 
that their X-ray emission is predominantly formed in the op¬ 
tically thick media. Similar to black holes, at lower luminosi¬ 
ties, log(Lx) ^ 36.5 — 37, n eutron stars u ndergo a transition 
to the hard spectral state fe.g. lBarreliZOOlll . The energy spec¬ 
tra in this state point at the low optical depth in the emission 
region, indicating a significant change of the geometry of the 
accretion flow. 

In the soft spectral state, the commonly accepted pic¬ 
ture of accretion at not too extreme values of accretion rate 
has two main ingredients - the accretion disk (AD) and the 
boundary layer (BL). While in the disk matter rotates with 


nearly Keplerian velocities, in the boundary layer it decel¬ 
erates down to the spin frequency of the neutron star and 
settles onto its surface. For the typical neutron star spin 
frequency 500 — 700Hz) compar able amounts of en- 

S r re released in these two regio ns dSunvaev & Shaku^ 
ISibgatullin & SunvaevI IZOOOl) . This picture is based 
on rather obvious qualitative expectations as well as 
more sophistic ated theoretical considerat i ons and numeri¬ 
cal modeling jSunvaey & Shakural Il986t iKluzniakl Il988 ' 


llnogamov & SunvaevI 1999t ISibgamllin & SunvaevI IzOOol) . 


It has been receiving, however, little direct observational con¬ 
firmation. Due to similarity of the spectra of the accretion 
disk and boundary layer the total spectrum has a smooth 
curved shape, which is difficult to decompo s e into sepa¬ 


rate spectral components \ Mitsudaeta illl984IWhite et al.L 
ll988tlDi Salvo et akLIZOOl : lDone et ^ 20021) . This made ap¬ 


plication of physically motivated spectral models to the de¬ 
scription of observed spectra of luminous neutron stars dif¬ 
ficult, in spite of very significant increase in the sensitivity 



























2 FOURIER-FREQUENCY RESOLVED SPECTROSCOPY OF LUMINOUS LMXBS 


2.1 The method 





Frequency. Hz 


Fig. 1. Average and frequency resolved spectra of GX340+0 
on the horizontal branch of the color-color diagram. The 
solid lines show the Comptonization spectrum with param¬ 
eters similar to those given in the Section^] 


of X-ray instruments. Not surprisingly, the best fit parame¬ 
ters derived from the data of different instruments and, corre¬ 
spondingly, the inferred values of the physically meaningful 
quantities are often in contradiction to each other. 

This ambiguity can be resolved if the spectral in¬ 
formation i s analysed to g ether with t iming data. Early 
results of iMitsiida et al.l (fl 984li and iMitsiida & Tanakal 
(Il986l) suggested that the boundary layer and accre¬ 
tion disk may have different patterns of spectral vari¬ 
ability. Based on the TENMA data, they studied the 
difference between the spectra averaged at different in¬ 
tensity levels - that rest ricted the range of accessible 
time scales to ^ 10^ sec. iGilfanov, Revnivtsev & Molkov 


(l2003l) and iRevnivtsev & Gilfanovri(2005l) have exploited 
the technique of Fouri er frequency resolved spectroscopy 
dRevnivtsev etallll999l) to study spectral variability of lumi¬ 
nous LMXBs in a broad range of time scales, including kHz 
QPO. Their findings are reviewed and discussed below. 


2. Fourier-frequency resolved spectroscopy of 
luminous LMXBs 


Fig. 2. Phase lags in GX340H-0 on the horizontal branch 
as function of energy {upper panel) and Fourier frequency 
{lower panel). The energy dependent phase lags were com¬ 
puted in the 1-32 Hz frequency range, the frequency de¬ 
pendent lags are between 3-6.5 keV and 6.5-13 keV energy 
bands. The phase is normalized to 0-1 interval. 


conventional (i.e. response folded) spectral approximations 
in order to describe the energy dependence of aperiodic vari¬ 
ability. Although the interpretation of the frequency resolved 
spectra is not always straightforward, several applications of 
this technique to var iability of black hole bin a ries gave mean¬ 
ingfu l results (e.g. IRevnivtsev et all 1 19991 IGilfanov et al.l 
120001 ) . 

We will use the following example as an illustration. 
Fet us consider a two-component spectrum, consisting of a 
constant and variable component. The variable component 
changes its normalization but not the spectral shape. In this 
case the shape of the frequency resolved spectrum would not 
depend on the Fourier frequency and would be identical to the 
spectrum of the variable component. Importantly, the X-ray 
flux in all energy channels will vary coherently and with zero 
time/phase lag between different energies. Presence of signif¬ 
icant phase lag and/or Fourier frequency dependence of the 
frequency resolved spectra would indicate that a more com¬ 
plex pattern of spectral variability is taking place. 


2.1. The method 

As defined in iRevnivtsev et all (Il999l) . the Fourier frequency 
resolved spectrum is the energy dependent rms amplitude 
in a selected frequency range, expressed in absolute (as op- 
posite to fract i onal) units. A similar approach was used by 
iMendez et alJ ll 19971) to study the energy spectrum of kHz 
oscillations in 4U0614-)-09. One of its advantages over frac¬ 
tional rms-vs.-energy dependence is the possibility to use 


With few exceptions ilDieters et all 1200(1) . phase lag be¬ 
tween light curves in different energy bands in luminous 
FMXBs is usually small, A(^ few x 10~coherence 
is consistent with unity, (e.g. IVaughan et al.[ 1 1994 Il999t 
IPieters et all 1200(1) and the behavior of the fractional rms- 
vs-e nergy dependence is simila r at different Fourier frequen¬ 
cies jvan der Khi Il986l I 2 OOOI) . This suggests, that Fourier 
frequency resolved spectra can be interpreted in a straightfor¬ 
ward and model-independent manner. 


































































3 BOUNDARY LAYER AND ACCRETION DISK EMISSION 


2 . 2 . Results 


For the case study we use archival data of PCA observations 
of a Z-source GX340-0 on the horizontal branch of the color- 
color diagram. Gilfanov. Revnivtsev & Molkov (120031) con¬ 
ducted similar study of an atoll source 4U1608-52 and arrived 
at similar conclusions. Our choice was defined by the require¬ 
ment that the PCA configuration combined sufficient energy 
resolution (large number of the energy channels) with good 
timing resolution and large total exposure time. The Fourier 
frequency resolved spectra in several frequency bands cor¬ 
responding to the band limited continuum noise component 
and the ~ 25 Hz QPO are shown in Fig. ^ along with the 
conventional spectrum of the source averaged over the same 
data. The figure clearly demonstrates that shape of the spectra 
depends on the Fourier frequency at low frequencies and be¬ 
comes independent of the frequency at / ^ 0.5 Hz. Another 
conclusion from the data presented in Fig. [2 important for 
the following discussion, is that all frequency resolved spec¬ 
tra are significantly harder than the average source spectrum. 

The phase lags as function of energy and Fourier fre¬ 
quency are shown in Fig. |2] No statistically significant phase 
lags were detected with an upper limit of A(p ~ 10“^, where 
phase (j) is normalized to the interval 0-1 (as opposed to 
0 - 27r). 


2.3. Interpretation 

We show below that independence of the frequency resolved 
spectra on the Fourier frequency and the smallness of the 
phase lags require a particularly simple form of the spectral 
variability. 

The constancy of the spectral shape with Fourier fre¬ 
quency implies that the power spectrum P{E, uj) can be rep¬ 
resented as a product of two functions, one of which depends 
on the energy and the other on the frequency only. For con¬ 
venience we write P{E, uj) in the form: 

P{E,u:)=S\E)x f{u:) (1) 

where non-negative functions S{E) and /(w) can be directly 
determined from the frequency resolved spectra. The Fourier 
image of the light curve E{E, t) is: 

F{E,u;) = SiE) X /(w) x (2) 

In the general case the complex argument (j){E, oj) can de¬ 
pend both on Fourier frequency w and energy E. If the phase 
lags between different energies are negligibly small, (j) de¬ 
pends on the Fourier frequency only and the Fourier image of 

F{E, t) is: 

F{E,uj) = S{E) X /(w) X (3) 

The light curve F{E, t) can be computed via inverse Fourier 
transform of F{E, uj): 

F{E,t) = J (EjF[E,uj)P‘^^ = 

= S{E) X [ = (4) 



Fig. 3. GX340-I-0: Coherence between the light curves in the 
3-6.5 and 6.5-13 keV energy bands as function of frequency. 
No correction for the dead time effects has been made. 


An arbitrary function of energy can obviously be added to the 
above expression: 

F{E,t) = So{E) + f{t)xS{E) (5) 

Thus, the light curves at different energies are related by a 
linear transformation. 

From Eq. it follows that the coherence of the signals in 
any two energy bands is exactly unity, as they are related by a 
linear transformation. This prediction is in a good agreement 
with observations (Fig.|3- 


3. Boundary layer and accretion disk emission 


As follows from eq. (Ell, two components can be distinguished 
in the GX340-0 emission. The term So{E) in Eq. (|3 is the 
constant (non-variable) part of the source emission spectrum 
and f{t) represents flux variations of the second, variable 
component.' The spectrum S (E) of the variable component 
does not change in the course of flux variations and equals the 
frequency resolved spectrum, i.e. can be directly determined 
from observations. 

There are two major components of accretion onto 
a slowly rotating weakly magnetized neutron star - 
(i) the Keplerian accretion disk and (ii) the bound¬ 
ary or spreading layer near the surface of the neu¬ 
tron star, in which the accreting matter decelerates 
to the spin frequency of the star a nd spread s over 
its surface (ISurivaev &^hakural 122^1_ Kjiizniakl Il22^ 
llnogamov & Siinvaevl ll QQol Popham feSunvaeu 120011) . 
These two geometrically distinct regions give com- 
parable contributions to the observed X-ray emission 
ISibgatullin & Sunvaevl l2000l) . Recalling eq.)|3, it is plau¬ 
sible to assume, that the variable part of the X-ray emis¬ 
sion is associated with one of these components. In order to 
check this assumption and to identify the variable component 
we consider below theoretical expectations for the disk and 


* The possibility of small variations of the spectral pa- 
rameter (e.g. temperature, opt i cal d epth etc.) is discussed by 
iGilfanov. Revnivtsev & MolkovI i2003h and is shown to be incon¬ 


S{E) X fit) 


sistent with observations. 


























4 BOUNDARY LAYER SPECTRUM 



Fig. 4. The average and frequency resolved spectra of 
GX340+0 (horizontal branch). The shaded area shows 
the plausible range of the boundary layer spectra calcu¬ 
lated as described in the section |3 The dashed (blue) 
histogram shows the best fit accret ion disk spectrum 
dOilfanov, Revnivtsev & Molkovl 1200.31) . The upper solid 
(red) histogram shows the boundary layer spectrum com¬ 
puted as the difference between the (observed) total and (pre¬ 
dicted) accretion disk spectrum. The lower solid histogram is 
the same but scaled to the total energy flux of the frequency 
resolved spectrum. 


boundary layer spectra and compare them with the observed 
frequency resolved spectra. 


At sufficiently high values of M both BL and disk are 
optically thick, as confirmed by the softness of the LMXB 
spectra. Simple arguments, taking into account the difference 
in the emitting areas suggest that the spectrum of the bound¬ 
ary layer should be harder than that of the accretion disk 
fe.g. lMitsuda et allll984tl&ebenev & Sunvaevll2002l) . Due 
to complexity of the boundary/spreading layer problem the 
theory has not advanced significantly beyond this qualitative 
statement - no models capable to directly predict its spec¬ 
trum exist yet. Significantly better progr ess has been achieved 
in modeling spectra of accretion di s ks llShakura & Sunvaevl 
I 1973 I: IShimura & Takaharal Il995t IRoss & Fabian! I 1994 . 


Relatively simple models of multicolor disk type which ac¬ 
count for the effects of Compton scattering with a simple 
color-to-effective temperature ratio turned out to be success¬ 
ful in describing the accretion disk sp ectra observed in the 
high state of black hole systems (e.g. lEbisawa et all Il991t 
lGierlin^ietal.[rT^ . 


For this reason we chose to use a model of the disk emis¬ 
sion as the starting point. The BL spectrum is computed as 
a difference between the (observed) total spectrum and the 
(predicted) disk spectrum. To estimate the plausible range of 
the BL spectra we investigate the parameter space of the ac¬ 


cretion disk model. For the latter we adopt the genera l rel¬ 
ativistic accretion disk model bv iLbisawa et alJ i(l99ll) (the 
“grad model” in XSPEC). The parameters of the model are; 
the source distance D, mass of the central object Mns, disk 
inclination angle i, the mass accretion rate M and the color- 
to-effective temperature ratio / = Tcoi/Tes- With this ap¬ 
proach we can predict the disk and the boundary layer spec¬ 
tra based on the observed X-ray flux and spectrum and very 
generic system parameters, such as neutron star spin fre¬ 
quency, the source distance etc. The procedure is described 
in detail in Gilfanov, Revnivtsev & Molkov il20n3h . The ob¬ 
tained range of the BL spectra is shown in Fig. |3 as the 
shaded area. The similarity of the predicted BL spectrum and 
of the observed frequency resolved spectrum is obvious. On 
the other hand the disk spectrum is significantly softer and is 
inc onsistent with the frequency reso l ved sp ectrum. 

Gilfanov. Revnivtsev & Molkovi il2flfl3l) conducted simi¬ 
lar investigation of an atoll source 4U1608-52 having 10 
times lower luminosity. They considered frequency resolved 
spectra at different frequencies, including the kHz QPOs and 
showed that it demonstrates behavior identical to GX340-0. 

Based on these results we conclude that the X-ray vari¬ 
ability on ~ second - millisecond time scales is related to 
variations of the luminosity of the boundary layer. The shape 
of its spectrum remains nearly constant in the course of these 
variations and equals the frequency resolved spectrum, i.e. 
can be directly obtained from the observations. This can be 
used to separate the boundary layer and the accretion disk 
contribution to the total spectrum and permits to check quan¬ 
titatively the predictions the accretion disk and boundary 
layer models. It also opens the possibility to measure rela¬ 
tive contributions of these two components of the accretion 
flow to the total observed X-ray emission. 


4. Boundary layer spectrum 

Based on the assumption that the frequency resolved spec¬ 
tra in luminous LMXBs adequately represent the boundary 
layer spectrum, we compare several LMXBs - atoll sources 
4U1608-52 and 4U1820-30 and Z-sources Cyg X-2 and GX 
17 h- 2. Their frequency resolved spectra at frequencies / ^ 
few Hz are shown in Fig|3 As before, for Z sources we used 
only data on the horizontal branch of the color-color dia¬ 
gram, where the amplitude of variability at these frequencies 
is maximal. To facilitate comparison, the normalizations of 
all spectra were adjusted to match that of GX340H-0. Sim¬ 
ilarity of the spectra is remarkable, especially considering 
significant difference in the average spectra and a factor of 
~ 10 — 20 spread in the luminosity between atoll and Z- 
sources (^ O.lMEdd and ^ MEdd correspondingly). 

The independence of the spectrum of the boundary layer 
on the luminosity lends support to the theoretical predictions 
bv rinogamov & Siinvaevlil 19991) that the boundary layer is ra¬ 
diation pressure supported, i.e. radiates at the local Eddington 
flux limit. In this model the luminosity of the spreading layer 
on the surface of the neutron star changes due to variations 
of its area, rather than surface emissivity of the unit area. If 
this picture is correct, the parameters of the BL emission can 









































5 NATURE OF THE Z-TRACK 


be used to determine the value of the Eddington flux limit on 
the surface of the neutron star. As the Eddington flux limit 
is uniquely determined by the neutron star surface gravity 
and the atmospheric chemical comp osition, the neutron star 
mass and radius can be constrained ( IRevnivtsev & Gilfanovl 
l2005h . 

The similarity of the spectral shape of the BE spectrum in 
different sources (Eig|3 indicates that there is no significant 
spread in the values of the mass and radius among LMXBs, 
in particular, that the surface gravity in atoll and Z-sources is 
similar. It also shows that there are no significant differences 
caused by variations in the atmospheric chemical abundances 
between sources. In particular we did not And statistically 
significant difference between ultra-compact compact binary 
4U1820-30 and other sources. 

The shape of the frequency resolved (ss boundary layer) 
spectrum can be adequately described by the saturated Comp- 
tonization. Eor the sake of comparison with other results and 
for convenient parameterization of the BL spectrum we used 
the Comptonization model of Titarchuk (1994). The best fit 
parameters of the model fitted in the 3-20 keV range simulta¬ 
neously to all five spectra shown in the Eig|3are: temperature 
of seed photons kTg = 1.5 ± 0.1, temperature of electrons 
kTe = 3.3 ± 0.4 and the optical depth r = 5 ± 1 for slab 
geometry. The best fit model is shown by the thick dotted 
line on Eig|3 The temperature of the black body spectrum 
describing the shape of the cutoff in the observed spectrum at 
energies >13 keV is /cTbb = 2.4 ± 0.1 keV (thin dashed line 
on Eig|3- 

The fact that kHz QPO show the same behavior as other 
components of the aperiodic variability indicates, that they 
have the same origin, i.e. are caused by the variations of the 
luminosity of the boundary layer. Although the kHz “clock” 
can be in the disk or due to it’s interaction with the neutron 
star, the actual modulation of the X-ray flux occurs on the 
neutron star surface. 

4.1. BL spectrum on the normal branch 

Eurther along the Z-track of GX340H-0, on normal and flar¬ 
ing branches, the fractional rms of the X-ray variability de¬ 
creases significantly, by a factor of ~ 5 — 10. Nevertheless, 
the statistics is sufficient to place meaningful constrains on 
the first half of the normal branch. The data indicates that the 
behavior of the frequency resolved spectra does not change 
its character - at sufficiently high frequency, / ^ 1 Hz, their 
shape does not depend upon the Eourier frequency and is 
significantly harder than the average spectrum and expected 
spectrum of the accretion disk. Therefore, we can conclude 
that frequency resolved spectra are representative of the spec¬ 
trum of the boundary layer. Eit to the frequency resolved 
spectrum by Comptonization model requires infinitely large 
values of the Comptonization parameter. Correspondingly, it 
can be described by Wien or blackbody spectrum (they are 
close to each other E ^ 3 keV range) with the best fit tem¬ 
perature of kT ss 2.4 keV. The frequency resolved spectra 
(saboundary layer spectra) of GX340-0 on the normal and 
horizontal branch are compared in Eig.|^ Thus, with increase 



Fig. 5. Eourier-frequency resolved spectra (saboundary 
layer spectra) of 5 Z- and atoll sources (from 
IRevnivtsev & Gilfanovl l2005l) . Eor 4U1608-52 the fre¬ 
quency resolved spectrum of the lower kHz QPO is shown. 
All spectra were corrected for the interstellar absorption. The 
thick short-dashed line shows the best fit Comptonization 
model with kTs = 1.5, kT^ = 3.3 keV, r = 5. The thin long- 
dashed line shows the blackbody spectrum with temperature 
fcTbb = 2.4 keV. 


of the mass accretion rate up to a value close to critical Ed¬ 
dington rate the boundary layer spectrum in the 3-20 keV 
energy range approaches a Wien spectrum. 

It would be interesting to follow up on these results and 
to consider the change of the BL spectrum from the hori¬ 
zontal to normal branch in other Z-sources. Unfortunately, in 
other four Z-sources from our sample the variability level on 
the normal branch is insufficient to obtain frequency resolved 
spectra with reasonable signal-to-noise ratio. 


5. Nature of the Z-track 


Knowledge of the shape of the boundary layer spectrum al¬ 
lows us to resolve the degeneracy caused by the similarity 
of the accretion disk and boundary layer spectra, which hin- 
dered many previous LMXB studies. As demonstrated by 
Gilfanov, Revnivtsev & Molkov ll2003ll . the spectra of atoll 
and Z-sources can be adequately described by the sum of 
the (renormalized) frequency resolved spectrum, represent¬ 
ing the boundary layer component, and of the accretion disk 
emission (Eig.l^. The spectrum of the latter is well described 
by the general relativistic accretion disk model. The best fit 
values of the mass accretion rate are consistent with those in¬ 
ferred from the observed X-ray flux and accretion efficiency 
appropriate for a 1 .4 Mq neutron star with spin frequency of 



































5 NATURE OF THE Z-TRACK 


5.1 BE fraction 
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Fig. 6. The absorption corrected frequency resolved spectra 
of QPO (~ boundary layer emission) in GX340 on the hori¬ 
zontal branch, (lower M) and upper half of the normal branch 
(higher M). The horizontal branch data is same as in Fig. [D- 
0] The solid line shows Wien spectrum with kT = 2.4 keV. 


Fig. 7. Dependence of the boundary layer contribution to the 
total X-ray emission of Z-sources as a function of the position 
on the Z-diagram 


~ 500 Hz. The agreement is especially remarkable, as the lu¬ 
minosity and mass accretion rate in atoll and Z-sources differ 
by the factor of ~ 10. 

We further exploit this approach and study the behavior 
of Z-sources in the color-color diagram in an attempt to re¬ 
late the motion along the Z-track to changes of the physi¬ 
cally meaningful parameters. We consider spectra integrated 
over 128-sec time intervals. As above, these spectra are fitted 
with a model, consisting of the boundary layer and the accre¬ 
tion disk components. The shape of the boundary layer spec¬ 
trum was approximated by the comptt model with param¬ 
eters from the section For the accretion disk spectrum we 
adopt the multicolor disk model (diskbb model in XSPEC) or 
general relativistic accretion disk model (grad). The model 
adequately describes observed spectra on the normal and hor¬ 
izontal branch. Further details of the analysis method, limita- 
tions of the model and discus sion of results are presented in 
iRevnivtsev & Gilfanovl(l2005l) . 

5.1. BL fraction 

The dependences of the BL contribution to the total X-ray 
emission on the position on the Z-track are plotted in Fig0 
The coordinate along the Z-track was defined to be propor¬ 
tional to the hard color with the reference points Sz = 1,2 
corresponding to the turning points of the “Z”, as is com¬ 
monly used for such plots. Statistical uncertainties in the 
values of the BL fraction are small and can be neglected, 
as confirmed by the dispersion of the points in Lig0 More 
important are the systematic ones associated with the im¬ 


precise knowledge of the shape of the BL spectrum and its 
possible variations along tl^ Z-track. As demonstrated in 
IRevnivtsev & GilfanovI i&f)^ their amplitude does not ex¬ 
ceed ~ 0.05 — 0.1 in the units of Lig.0 

The Lig0 suggests that the boundary layer fraction de¬ 
creases along the Z-track and it is smaller on the nor- 
mal branch than on the horiz ontal branch. As discussed in 
IRevnivtsev & GilfanovI ll2005l) . this conclusion is rather ro¬ 
bust, as long as the assumption regarding the constancy of 
the boundary layer spectrum is approximately correct. As 
the variability at / ^ 1 Hz is primarily associated with the 
boundary layer emission, the decrease of the boundary layer 
fraction along the Z-track also explains well-known decrease 
of the level of aperiodic and quasi-periodic variability. 

Although no simple physical interpretation of the ob¬ 
served behavior can be offered, we mention several possibili¬ 
ties. One of these is that the general structure of the accretion 
flow does not change significantly and ~ 50% of the energy 
is always released on, or very close to the neutron star sur¬ 
face. The apparent decrease of the boundary layer fraction 
on the normal branch is a result of its geometrical obscura¬ 
tion by, for example, the geometrically thickened accretion 
disk. An alternative possibility is that at high values of the 
mass accretion rate M ~ MEdd a significant modification of 
the accretion flow structure occurs and its division into two 
geometrically distinct parts - boundary layer and accretion 
disk, becomes inapplicable. Namely, due to non-negligible 
pressure effects the deceleration of the orbital motion of the 
accreting matter from Keplerian frequency to the neutron star 
spin frequency would take place in a geometrically extended 
region with the radial extend of AR ~ i?NS- In this case, the 

























6 SUMMARY 



Soft color (4-6 keV)/(3-4 keV) 


Fig. 8. The horizontal and normal branches in the color-color 
diagram of several Z-sources (Cyg X-2, GX340H-0, Sco X-1, 
GX 5-1 and GX 17 h- 2). Overlayed on the data are the tracks 
predicted by the model (section l5^ . The two thick dashed 
lines show evolution of the colors with change of the accre¬ 
tion rate in the disk for two different values of the BL frac¬ 
tion: 44% (upper) and zero (lower) The values of the mass 
accretion rates M of the disk component are marked in units 
of 10^® g/s. The thick solid line shows the Z track with tran¬ 
sition at M ~ 2 ■ 10^® g/s. 


observed decrease of the boundary layer fraction could reflect 
actual decrease of the fraction of the energy released on the 
neutron star surface with the rest of the energy being released 
in the extended transition region. 


5.2. Shape of Z-diagram 

Motivated by these results we use the two-component spec¬ 
tral model to explain the shape of the Z-track on the color- 
color diagram. For this purpose, the accretion disk compo¬ 
nent is modeled by the general relativistic accretion disk 
model of Ebisawa et al. (1991) (the grad model in XSPEC) 
which explicitly includes dependence on the mass accretion 
rate. Using this spectral model we can calculate the position 
in the color-color diagram as a function of the mass accretion 
rate and the boundary layer fraction. This is an attempt to un¬ 
derstand general tendencies of the Z-diagram, rather than to 
construct its precise quantitative description. 

The results are presented in Eig|3 The two dashed curves 
in the figure show the evolution of spectral colors with the 
increase of M for two values of the boundary layer fraction, 
BL/disk=0.8 and BL fraction of zero. The evolution of col¬ 
ors corresponding to a change of the BL contribution from 


TsL/Tdisk = 0.8 to zero at M ~ 2 x 10^® g/s is shown by 
the thick solid line. 

The general shape of the Z-track can be reproduced in the 
model as a result of variation of two parameters - the mass 
accretion rate and BL fraction. The mass accretion rate in¬ 
creases along the Z-track. The Z-shape of the track is defined 
by the variation of the BL fraction, which decreases along 
the normal branch from the value of ~ 50% expected in the 
“standard” theories to a small number of the order of ~zero 
at the end of the normal branch. The exact value of M, cor¬ 
responding to the transition from the horizontal to the normal 
branch depends on the disk model parameters - the binary 
system inclination, the mass of the neutron star and the spec¬ 
tral hardening factor. Lor our choice of parameters, it equals 
M ~ 2 ■ 10^® g/sec, i.e. is of the order of the Eddington crit¬ 
ical value for a 1.4M0 neutron star. 


6. Summary 


1. The X-ray variability in luminous LMXBs on the short 
timescales, / ^ 1 Hz, is caused by variations of the lumi¬ 
nosity of the boundary layer. The accretion disk emission 
is significantly less variable at these frequencies. The BL 
spectrum remains nearly constant in the course of lumi¬ 
nosity variations and its shape equals the frequency re¬ 
solved spectrum, i.e. can be directly derived from the tim¬ 
ing data (Ligl^l. 

2. In the investigated range of the mass accretion rate M ^ 
(0.1 — l)MEdd. the boundary layer spectrum depends 
weakly on M. Its shape is remarkably similar in atoll and 
Z-sources (Fig. |3i, despite an order of magnitude differ¬ 
ence in the mass accretion rate. Data indicates that in the 
limit of high M ~ MEdd. the boundary layer spectrum 
can be described by Wien spectrum with kT Ri 2.4 keV 
(Fig-E). At lower values of M the spectra are better de¬ 
scribed by model of saturated Comptonization with elec¬ 
tron temperature of ~ 2 — 4 keV and Comptonization 
parameter y ^ 1. Weak dependence of the BL spectrum 
on the global valu e of M lends support to the theoret¬ 
ical suggestion by llnogamov & ■Siinvaevl ill999l) that the 
boundary layer is radiation pressure supported. 

3. The kHz QPOs appear to have the same origin as ape¬ 
riodic and quasiperiodic variability at lower frequencies. 
The msec flux modulations originate on the surface of the 
neutron star although the kHz “clock” might reside in the 
disk or be determined by the disk - neutron star interac¬ 
tion. 

4. We attempt to relate the motion of Z-sources along the Z- 
track to changes in the values of the physically meaning¬ 
ful parameters. Our results suggest that the contribution 
of the boundary layer component to the observed emis¬ 
sion decreases along the Z-track from the conventional 
value of ~ 50% on the horizontal branch to a rather small 
number at the end of the normal branch (Fig l7l8> . The 
main difference of our approach from previous attempts 
is in the a priori knowledge of the shape of the bound¬ 
ary layer spectrum. This allowed us to avoid ambiguity of 
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the spectral decomposition into boundary layer and disk 
components. 
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